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Abstract: We present the NLO Susy-QCD corrections to the production of a pair of the 
lightest neutralinos plus one jet at the LHC, appearing as a monojet signature in combi- 
nation with missing energy. We fully include all non-resonant diagrams, i.e. we do not 
assume that production and decay factorise. We derive a parameter point based on the 
pl9MSSM which is compatible with current experimental bounds and show distributions 
based on missing transverse energy and jet observables. Our results are produced with the 
program GoSam [1] for automated one-loop calculations in combination with MadDipole/- 
MadGraph for the real radiation part. 
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1. Introduction 



With the LHC experiments performing extremely well, we have entered the exciting phase 
where we can investigate the properties of a new discovery [2, 3], and get exclusion bounds 
and hopefully also hints for Beyond the Standard Model physics at energy ranges which 
have never been probed before. 

If the new boson with a mass around 125 GeV is a fundamental scalar, we have to 
figure out why it is protected from higher order corrections of the order of the Planck scale. 
Supersymmetry (for a review see e.g. [4, 5] and references therein) offers an elegant expla- 
nation for a stabilisation mechanism, and in addition contains massive weakly interacting 
particles which can serve as dark matter candidates. 

While the hopes of an early Susy discovery at the LHC have withered as recent LHC 
measurements have pushed up the lower limits on squark and gluino masses considerably [6, 
7, 8, 9, 10, 11], there is no stringent lower limit on the mass of the lightest neutralino [12, 13], 
and the pair production of charginos/neutralinos becomes increasingly important as a 
"discovery channel" of Supersymmetry. 

In many Susy scenarios, the neutralino x? is assumed to be the lightest supersymmetric 
particle (LSP) and thus is stable if R-parity is conserved. Therefore, XiXi production, 
either directly or through the decay of heavier neutralinos/charginos, leading to signatures 
of missing energy and energetic jets/leptons, is a process of primary interest in the context 
of current Susy searches. Hence it is desirable to have predictions which include the NLO 
Susy-QCD corrections to such processes, not only at the level of total cross sections, but also 
for differential distributions. While neutralino pair production without any additional jets 
is not very illuminating from an experimental point of view, the production of neutralino 
pairs in association with additional jets is interesting, since the signature "missing energy 
+ jets" smells like New Physics, and the distributions of jet observables can be used to 
investigate the nature (e.g. the spin) of the object(s) carrying away the missing energy. 

A very clean signal of new physics would be the observation of an excess in events 
involving a very energetic monojet in combination with missing transverse energy. Searches 
for monojets at the LHC at 7 and 8 TeV have been carried out already [14, 15, 16, 17], 
and turned mainly into limits on models involving extra dimensions, assuming that the 
missing energy is due to graviton production. If the missing energy is due to neutralinos, 
studying such processes could provide information on the nature and couplings of the LSP 
and thus on dark matter [18, 19]. Further, monojet signatures are interesting in the context 
of constraints on invisible decays of the Higgs boson, as the invisibly decaying boson may 
recoil against hard QCD radiation [20, 21, 22]. 

Total cross sections for the production of charginos and neutralinos at next-to-leading 
order in QCD have been calculated in [23]. Recently, updates for the LHC at 7 and 8 TeV 
and current popular benchmark points [24] have been given in [25, 26]. Resummation of 
large logarithms in the threshold and small-pT regions has also been carried out [27, 28, 29] 
for the production of gaugino pairs. However, predictions for neutralino pair production 
in association with one or several jets in a differential form were, until recently, only avail- 
able at leading order [18]. The first differential NLO QCD corrections for the final state 



- 1 - 



of missing transverse energy in association with two jets, where the missing energy stems 
from the decay of a squark pair into quarks and neutralinos, have been calculated in [30]. 
For the case of squark and gluino pair production without decay or additional jets, the 
number of available results beyond the leading order is larger than for charginos/neutrali- 
nos, because these processes were hoped to be seen already at the Tevatron or at early 
stages of LHC measurements. The first NLO calculations [31, 32, 33, 34], partly enter- 
ing the code Prospino [35], were followed by electro- weak corrections [36, 37, 38, 39, 40], 
resummation [41, 42, 43, 44, 45] and NNLO threshold corrections [46, 47]. NLO QCD cor- 
rections to squark and gluino pair production compared to results from LO matrix element 
plus parton shower merging have been presented recently in [48], based on calculations 
in [49, 50]. 

Next-to- leading order predictions involving chargino/neutralino pairs used in experi- 
mental analysis were usually obtained by calculating the NLO K-factors for the total cross 
sections using e.g. the code Prospino [35, 33], and rescaling the LO predictions accord- 
ingly. However, it is not at all guaranteed that the same K-factors can be applied globally, 
in particular in the presence of stringent search cuts. QCD radiation can change the shape 
of the distributions and affect the mass and spin measurements considerably [51]. If the 
supersymmetric spectrum is highly compressed, QCD radiation can also seriously affect 
the exclusion bounds extracted on the basis of leading order assumptions [21, 52]. 

In this paper, we calculate the NLO QCD corrections to the production of a pair of the 
lightest neutralinos in association with one jet. We include not only resonant contributions 
from squark decay, but also all non-resonant contributions. We present our results in a 
fully differential form, showing distributions for observables involving the jet and missing 
energy, which can be compared straightforwardly to data. The treatment of diagrams 
involving resonant squarks needs special attention, as the NLO real corrections formally 
also contain diagrams which can be regarded as leading order contributions to a different 
process, which is resonant squark pair production and subsequent decay. We also calculate 
contributions from Higgs production through a heavy quark or squark loop. However, 
these contributions are found to be numerically very small. For our studies we consider the 
phenomenological MSSM (pMSSM) [24, 53, 54], in a variant involving 19 free parameters 
(pl9MSSM). In this framework, we derive a point where the mass of the lightest Higgs boson 
h is mh = 125.8 GeV. The virtual corrections have been calculated with the automated 
one- loop program GoSam [1], where the integrals involving complex masses have been 
called from the integral libraries Golem95 [55, 56] and OneLOop [57]. The real radiation 
matrix elements are generated using MadGraph [58] and MadDipole [59, 60]. 

The paper is organized as follows. In Section 2, we describe our method for the calcu- 
lation of the virtual and real corrections, in particular the generation of the renormalisation 
counterterms and the treatment of resonant squarks, the latter being further discussed in 
the Appendix. We also include a phenomenological discussion of MSSM parameter points. 
In Section 3 we present our numerical results, before we conclude in Section 4. 
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2. Calculational framework 
2.1 Virtual corrections 

The one-loop virtual contribution to the NLO result is calculated using the program Go- 
Sam [1]. We use FeynRules [61] to produce a model file in the UFO [62] format that can 
be read directly by GoSam. 

For the virtual amplitude we have 0(1400) diagrams to calculate for each subprocess. 
We neglect 6-quarks in the initial state. The most complicated diagrams are rank-3 pen- 
tagons, with up to 4 internal masses. We illustrate two of the pentagon diagrams in Fig. 1. 
We include finite widths in the loop integrals and therefore we need a basis set of complex 
integrals, which we call from the libraries Golem95C [56] and OneLOop [57]. To calculate the 
loop amplitude in a numerically robust way, GoSam is able to interchange between differ- 
ent reduction schemes at runtime. Our default reduction strategy is to use Samurai [63] 
and, if it fails, to reprocess the point with Golem95C [56] using tensorial reconstruction at 
integrand level [64] . Due to the large internal masses and the small squark widths present in 
the integrals, and due to the high rank of the pentagons, numerical stability is a nontrivial 
issue in this process. Therefore it is crucial that we have this rescue system available during 
the numerical integration. We use the dimensional reduction scheme (DRED) where only 
the internal momenta are kept in D dimensions. We also calculate contributions where 
neutral Higgs bosons can be produced by a loop-induced process. These can then decay to 
a pair of the lightest neutralinos. 




Figure 1: Two illustrative pentagon diagrams calculated by GoSam. 



2.2 Counter Term Diagrams 

Due to the appearance of self-energy insertions that can have internal masses different to 
the mass of the incoming particle we modified the handling of counter term diagrams from 
the default GoSam template file. In the default GoSam template file each self-energy 
insertion in the virtual diagrams is rewritten as a sum of the original part plus a part that 
integrates to give the counter term diagram. Here we pursue a different approach as follows: 
each diagram containing a self-energy insertion is "tagged" when the diagram topologies 
are analyzed. This information is passed through the processing of the diagrams and is 
available when the numerators of the diagrams are reduced. The counter term diagram 
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is then reconstructed from the original diagram, by shrinking the self-energy insertions 
to a point. This is then reduced alongside the loop diagram and written to a separate 
file for numerical evaluation. The counter term diagrams depend on model dependent 
renormalisation constants that we calculate separately. It is planned that this approach 
will be incorporated into a future release of GoSam. 

2.3 Real corrections and numerical integration 

The leading order (LO) and the real radiation matrix elements are generated using MadGraph 
[58]. For the subtraction of the infrared singularities we use Catani-Seymour dipoles [65], 
supplemented with a phase space restriction parameter a as proposed in [66] and imple- 
mented in the package MadDipole [59, 60]. The generation of the various pieces of the 
code (tree-level, real emission process, subtraction terms) and their combination with the 
integration routines was performed in a fully automated way. For the numerical integra- 
tion we used MadEvent [67, 68], slightly modified in order to integrate the different NLO 
contributions. 

2.4 Treatment of diagrams with resonant squarks 

For processes involving unstable particles, the proper definition of the set of diagrams 
contributing to the next-to-leading order corrections is not obvious. There are problems 
of double counting as diagrams with additional real radiation from the unstable particle 
in the final state can, if it becomes resonant, also be regarded as part of a leading order 
process with the decay already included in the narrow width approximation. The problem 
is similar to the case of Wt and W + W~ b production at NLO, where the W + W~ bb final 
state occurring in the NLO real corrections, if stemming from doubly resonant top decays, 
can also be viewed as belonging to leading order ti production and decay. This problem 
has been discussed in detail in [69, 70] for the case of Wt production. The case at hand is 
very similar, with t — > W b replaced by q — > x± Q- Consider for example the leading order 
diagrams in Fig. 2. Squark exchange in the t-channel, as shown in Fig. 2(a), cannot lead to 
any resonance, but in the case of s-channel squark exchange shown in Fig. 2(b), the squark 
can become resonant, and it can be viewed as a diagram for squark production in association 
with a neutralino, with squark decay included in the narrow width approximation. Now at 
NLO, when the real radiation of an additional parton is included, a new channel opens up, 
where two squarks can decay resonantly into a quark and a neutralino, as shown in Fig. 3. 
Close to the resonance, this contribution gets quite large, and in fact should rather be 
counted as a leading order contribution to squark pair production with subsequent squark 
decay, because here we are interested in the radiative corrections to the final state of a 
monojet in association with a neutralino pair. 

For this reason the calculation was carried out in two different ways. In the first 
approach we take into account all possible diagrams leading to the required final state con- 
sisting of two neutralinos and two QCD partons. In particular this includes the possibility 
of having two on-shell squarks. 

As it turns out, the real radiation part is vastly dominated by these diagrams, spoiling 
the convergence of the perturbation series. The result behaves like a tree-level calculation, 
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Figure 2: Examples of a leading order diagrams with a squark in (a) the t-channel, (b) the 
s-channel. 



involving one more order in the strong coupling constant than the tree level for the monojet 
case. The difference to an actual tree-level calculation is given by the fact that one parton 
can become unresolved, leading to the signal "single jet plus missing energy", while the 
infrared singularities due to the unresolved parton are canceled by the virtual corrections 
or absorbed into the parton distribution functions. Therefore, even though the calculation 
which includes these resonant diagrams cannot be regarded as a genuine next-to-leading 
order correction to the cross section describing a neutralino pair in association with one 
jet, it is still a physically meaningful quantity. However, from an experimental point of 
view, a complete description of the final state of missing energy plus up to two jets would 
be more useful. This however would require the full NLO calculation of the production of 
a neutralino pair in association with two jets, where the jets can either originate from the 
decay of squarks and gluinos or be produced directly from partons in the hard interaction. 
This is a very complex task which is beyond the scope of this paper. 

In the second approach we follow a strategy proposed in [69], namely we remove the 
diagrams with two squarks in the s-channel from the amplitude. Removing diagrams from 
the amplitude generally violates gauge invariance. Ref. [69] contains a study about the 
impact of violating gauge invariance by such a removal of diagrams, where the effects were 
found to be small for commonly used gauges. 

For our calculation we assume that the largest contribution of diagrams with two 
squarks in the s-channel come from those points in phase space where both squarks are on- 
shell, and that off-shell effects are suppressed by a factor of T/M. Therefore it is sufficient 
to consider the 2—7-2 process of producing two squarks. In this case one can show that 
the gauge dependence vanishes for covariant gauges and for a large class of non-covariant 
gauges. We give a proof in appendix A. 

Another solution to the double counting problem based on the subtraction of the res- 
onant diagrams is the so-called "prospino scheme" [33, 50], where in its most recent vari- 
ant [48] a subtraction term is introduced which involves a momentum remapping oriented 
at Catani-Seymour mappings to preserve the on-shell conditions for both intermediate and 
final state momenta. The subtraction term is then constructed in such a way that it comes 
into operation when the intermediate particle goes on-shell. However, for multi-particle 
final states where the Born process is already 2 -4 3 or of higher multiplicity, and doubly 
resonant diagrams appear at NLO, the generalisation of this procedure is not obvious. 
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Figure 3: Example of a "doubly resonant" squark diagram occurring in the NLO real corrections. 



2.5 Phenomenological discussion of Susy parameter points 

The discovery of a boson with mass of about 125 GeV [2, 3] and a lack of any experimental 
evidence for supersymmetry has put many Susy models under strain. In particular, the 
family of the so-called constrained MSSM (cMSSM) is increasingly unfavoured in light 
of the present experimental data [71, 72]. Therefore we choose a more pragmatic and 
experimentally motivated parameterisation of Susy, known as the phenomenological MSSM 
(pMSSM) [53, 54, 24], in a variant involving 19 free parameters (pl9MSSM). In the pMSSM, 
no underlying Susy breaking scenario is chosen, and gauge coupling unification is not 
postulated. However, it is assumed that there are no sources of CP violation and generation 
mixing beyond the Standard Model ones, so that all off-diagonal elements in the sfermion 
mass matrices are equal to zero, and the first and second-generation soft terms are equal. 
We will further assume that the LSP is the lightest neutralino. 

In such models we can corner the parameter space where the squarks, gluino and 
neutralinos have masses close together. This region is of particular interest for the process 
considered here, as the squark masses can be low enough for the signal to be significant at 
LHC energies while evading the current LHC exclusion regions. 

Compressed spectra of this type were also studied in [21, 52, 12] where the search for 
Susy through the recoil of light Susy particles against initial state radiation is discussed. 
Our choice of Susy parameter point can be viewed as qualitatively like the "Equal Mass" 
scenario in Ref. [52]. In this situation we expect to see soft jets from the decays of the 
squarks to the neutralino. 

For our results we choose a modification of the parameter point pl9MSSMl.l which we 
call the pl9MSSMlAmod. The pl9MSSMl line was introduced in [24] and was constructed 
as a benchmark for these compressed Susy scenarios, and is labelled by an integer N, 
i.e. pl9MSSMl.iV. The pl9MSSMl line is controlled by one parameter: the gaugino 
mass parameter Mi . The two lightest generation sfermion masses M £ 2 and the mass 
of the gluino, M 3 , are fixed at 1.2M X . pl9MSSMl.JV is defined for N = 1 at the value 
Mi = 300 GeV and this is increased by 100 GeV for each subsequent value of N. We 
effectively decouple all other particles in the model by setting the other mass parameters 
to a higher scale, here 2500 GeV. We choose a value of 10 for tan/3. 

Our choice modifies pl9MSSMl.l in two ways. Firstly, we make contact with the 
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point pl9MSSMlA, given in [30], by setting the heavier Higgs bosons mass inputs, \i and 
mAi to the higher scale. Like the heavier squarks these particles are effectively decoupled. 
Secondly, we further modify this point by choosing At = 5000 such that the mass of the 
lightest Higgs h in our model can be identified with the boson observed at the LHC with 
mass rrih = 125.8 GeV. We illustrate the effect that varying At has on the mass of the 
lightest Higgs boson from At = to the maximum Higgs mass in Fig. 4 over a range of the 
gaugino mass parameter M\ . We do not plot the theory uncertainty coming from unknown 
higher order corrections nor the uncertainty from the input parameters. The Susy masses 
that we use for our calculation, at this point in parameter space, are given explicitly in 
Table 2. 




105 

500 1000 1500 2000 2500 

< m- > [GeV] 

Figure 4: The lightest Higgs boson in the MSSM depends on the parameter At through one- 
loop and two-loop corrections. Here, the band illustrates how mt varies with At in the benchmark 
scenario pl9MSSMl. This variation is examined over a range of input parameters M\ which changes 
the spectrum of the light squark masses. The x-axis is the average value of the light squark masses. 

The particle spectrum was produced using SOFTSUSY [73] and then the decay widths 
calculated using the package SUSYHIT [74] in which we used the packages HDECAY [75] 
and SDECAY [76]. GoSam can read input cards provided in the SLHA format [77, 78] 
allowing one to change the Susy parameter point without recompiling the code. 

2.6 Higgs Contribution to the Signal 

As mentioned in Section 2.1 we also calculate the contribution arising from Higgs produc- 
tion through heavy quark and squark loops with subsequent Higgs decay to a neutralino 
pair. We illustrate these types of diagrams in Fig. 5. To quantify their contribution to the 
total cross section we can easily isolate these types of diagrams using the diagram filtering 
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system in GoSam. We find the Higgs boson contribution to the total cross section to be 
negligible, so these diagrams are not included in the results shown in Section 3. 




Figure 5: Examples of a loop diagrams involving the MSSM Higgs bosons. Here h stands for any 
of the MSSM Higgs bosons in the set {h, H, A}. 



2.7 Checks on the Result 



We have checked that after UV renormalisation, all poles from the virtual contributions 
cancel with the poles from the infrared insertion operator [65] in the real radiation. We have 
also checked the unrenormalised virtual matrix element against the program FeynArts/ 
FormCalc [79, 80, 81]. We found agreement for the partonic subprocess ug — > XiXi u - All 
other subprocesses can be found by exploiting crossing symmetry. 
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Figure 6: Cancellation of the finite remainders of the UV divergent parts of the amplitude with 
the corresponding counter terms in the resonance region. 



Furthermore, we checked our renormalisation procedure by investigating the cancella- 
tion of the finite remainders of the UV divergent loop contributions with the counter terms, 
as illustrated in Fig. 6. 
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3. Phenomenological results 

In this section we present a selection of phenomenological results for proton proton collisions 
at the LHC at 8 TeV. 

3.1 Setup and cuts 

For all the results and distributions shown in this section we have used the parameters 
listed in Tables 1 and 2. 



Standard Model Parameters 


M z = 91.188 GeV 
M t = 173.4 GeV 
a{M z ) = 1/127.934 


T z = 2.4952 GeV 
r t = 1.4384 GeV 



Table 1: Standard Model parameters used for the phenomenological study. 



SUSY Parameters 



M^o = 299.5 
M s = 415.9 
Mu L = 339.8 
M aR = 396.1 
M 3l = 348.3 
M 3r = 392.5 
M- H = 2518.0 
Mi r = 2541.8 

= 2403 • 7 
Mr = 2668.6 





= 


r§ = 


= 4.801 




= 0.002562 


Y- 


= 0.1696 




= 0.003556 




= 0.04004 


T l, 


= 158.1 




= 161.0 




= 148.5 




= 182.9 



Table 2: Masses and widths of the supersymmetric particles for the benchmark point used. The 
second generation of squarks is degenerate with the first generation of squarks. All parameters are 
given in GeV. 

The weak mixing angle is calculated from the W and Z masses. The strong coupling 
constant and its running are determined by the set of parton distribution functions. We 
used an NLO pdf set from NNPDF2.3 [82], where the values for a s at leading order and 
next-todeading order are given by 

a s {M z ) = 0.119 , 

and the running is calculated at one loop for the treedevel result and at two loops for the 
next-todeading order parts. As we neglect initial state 6-quarks, we use the Nf = 4 version 
of the pdf set. Further, we assume flavour-diagonal Susy-QCD couplings. 

For the jet clustering we used an anti-fer algorithm [83] with a cone size of R = 0.4 
provided by the Fast Jet package [84, 85]. We choose [i = Ht/2 for our central scale, where 
we define Ht as Ht = Yli ^T,i with i running over the momenta of the two neutralinos 
and the jet(s). 
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We use the following set of cuts 

p T (leadingjet) > 100 GeV, |r/j| < 4.5 . (3.1) 
In addition we impose a cut on the missing transverse energy of 

-St, miss 

> 85 GeV. (3.2) 

We also impose a jet veto of 30 GeV on a second jet which at NLO originates from the 
2 — > 4 part of the real radiation corrections, as discussed in Section 2.4. Our relatively 
low cut on the transverse missing energy is motivated by the fact that the neutralinos Xi 
could be rather light, and therefore the requirement of very large £/r,miss> which is well 
motivated in searches for graviton production in association with monojets, could be too 
restrictive in the case of neutralino pair production in association with one jet. In this 
case the neutralinos do not originate from long cascades of heavier objects with additional 
missing energy produced along the cascade. 

3.2 Numerical results 

In this section we show distributions for the observables p 3 T , the transverse momentum 
of the jet, the missing transverse energy E™ 1SS , and the angle (p(p miss , fp ') , where p mms 
is defined as minus the vector sum of the visible particles in the event, and j? is the 
momentum of the leading jet. We show two types of distributions for each observable: one 
where each distribution is normalized to one in order to exhibit the difference in shape, 
and another with absolute values. For the normalized histograms in Figs. 7, 9 and 11, we 
show results for both approaches, the one including the doubly resonant diagrams and the 
one with diagram removal. 

Investigating the behaviour of the cross sections under scale variations, we observe 
the following. As expected, the results including the doubly resonant diagrams show no 
improvement of the scale uncertainty at NLO, because they are completely dominated by 
the 2—7-4 real radiation, and therefore the scale dependence is not compensated by the 
virtual contributions. The case where these diagrams are removed is still dominated by the 
new channels opening up in the NLO real radiation contributions. Therefore, in this case 
we do not find a stabilisation of the scale dependence either. However, the cross sections 
are sizeable. Using fi = = fip and varying between Ht/^ < fJ> < 2Ht, we find 20 - 
30 fb for the LO cross section, while the NLO subtracted cross section amounts to about 
100 fb for the central scale, and the one including doubly resonant diagrams to 960 fb for 
the central scale. This means that the point pl9MSSMlAmod considered here could in 
principle be tested with the data accumulated so far. 

For the results including the doubly resonant diagrams, it is pointless to determine a 
K-factor, as in this case it is not well defined to which leading order process the higher 
order terms should be attributed, as explained in Section 2.4. 

For the subtracted case, where we show absolute as well as normalized results, the K- 
factors are still quite large, as can be seen from Figs. 8, 10 and 12. This can be understood 
as being mainly due to new channels opening up in the NLO real radiation part, in the 
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Figure 7: Normalized distributions showing the transverse momentum distribution of the leading 
jet for the process pp — > XiXi+j e t at y/s = 8 TeV, comparing the cases where the resonant diagrams 
are included to the ones where they are subtracted. 
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Figure 8: Transverse momentum distribution of the leading jet for the process pp — > XiXi+jet a t 
v/s = 8 TeV. 

presence of an additional QCD parton. The distributions also show that the K-factors are 
not uniform, which implies that the "standard" procedure to use K-factors based on total 
cross sections calculated at NLO can be misleading. For the angle 0(p miss , jP) between 
the jet and the missing momentum, shown in Figs. 11 and 12, NLO is obviously the first 
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Figure 9: Normalized distributions showing the missing transverse energy E™ lss , comparing the 
cases where the resonant diagrams are included to the ones where they are subtracted. 




order where the fixed order prediction makes sense, because at LO the vector sum of the 
neutralino momenta will always be back-to-back to the recoiling jet. 

The large K-factors can partly be attributed to the fact that at NLO, a considerable 
number of new partonic channels are opening up: the gluon-gluon initiated processes, and 
the ones initiated by uu, dd, ud, ud are all absent at leading order (where only qq and qg/qg 
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Figure 11: Normalized distributions showing the angle 4>(p lnlss ,jP), comparing the cases where 
the resonant diagrams are included to the ones where they are subtracted. 




Figure 12: The angle 4>(p inlss ,p') between the leading jet and missing momentum vectors, where 
the latter is defined as the negation of the visible momenta in the event. 

are present) and the new channels together make up almost 50% of the cross section. If 
we roughly estimate the K-factor which would result from partonic initial states which are 
already present at LO, it would amount to K ~ 2.3. Further, we investigated the point 
SPSla [86] just for reference, and found that for this point the K-factors are also smaller. 
This can be attributed to the fact that for the compressed spectrum we are considering here, 
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the widths of the first generation squarks are very small, leading to larger contributions if 
the squarks are close to being on-shell. 

4. Conclusions 

We have calculated the NLO Susy-QCD corrections to the production of a pair of the 
lightest neutralinos in association with one jet. We did not use the approximation of 
factorising production and decay, but fully included all non-resonant contributions. Con- 
tributions from Higgs production through a heavy quark or squark loop were calculated 
and found to be numerically negligible. The calculation has been performed using two 
different approaches to treat the doubly resonant diagrams appearing in the NLO real ra- 
diation contribution: one is based on diagram removal and the other based solely on a veto 
on the second jet. In the latter approach, the K-factors are obviously very large as the 
whole result is dominated by a reaction which can also can be viewed as the Born level for 
a different process (resonant squark pair production with a subsequent factorisable decay 
into a neutralino and a jet). We present our results in a fully differential form, based on 
the experimentally accessible jet and missing Et observables. 

On a technical level, to the best of our knowledge, this is the first 2—7-3 NLO cal- 
culation within the MSSM which includes full off-shell effects, as well as complex masses. 
Even though we only consider the phenomenological MSSM (pl9MSSM) here, we would 
like to emphasize that our setup is largely automated, using the public one-loop program 
GoSam in combination with MadGraph, MadDipole and FeynRules, such that other pa- 
rameter points, and even other models Beyond the Standard Model, can be studied as well 
within the same framework. 
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A. Gauge dependence 

In this appendix we examine the gauge dependence of the diagrams that have been removed 
from the amplitude in the real emission part as discussed in 2.4. We show that this gauge 
dependence vanishes for covariant gauges and for a large class of non-covariant gauges. 

The only diagrams, once omitted, that can lead to a dependence on the choice of 
gauge are of the type shown in Fig. 13. In this diagram there is an s-channel gluon which 
decays into a squark-antisquark pair. As the biggest contribution to the cross section 
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comes from the parts of the phase space where the two squarks are on-shell, it is sufficient 
for our argument to consider the 2—7-2 proccess of squark pair production and neglect the 
subsequent decay of the squarks. We denote the incoming momenta of the quarks as q\, q 2 
and the outgoing momenta of the squarks as p\,p 2 - in the following we neglect overall 
prefactors like color factors and coupling constants as they are irrelevant for our argument. 
The same holds for factors of i and any minus signs. The amplitude of the 2—)- 2 process 
can be written as 

M ~ v( qi ) 7 " IV u(q 2 ) ■ (tf - p») , (A.l) 
where D^ u denotes the gluon propagator, which in Feynman gauge is simply given by 



D 



' k 2 



with k = qi + q2 



(A.2) 



(A.3) 



(A.4) 



Choosing the Feynman gauge and contracting the Lorentz indices expression gives 

M ~ v(qi) (p/i -v/2) u(q 2 ) 
and after squaring and performing the fermion spin sum one obtains 

|A*| 2 ~tr(gi . 

As the gluon propagator is gauge dependent, 
the gauge dependence vanishes only in the sum 
of all contributing amplitudes. 

To calculate the effect of a specific gauge to 
the given diagram we start with a general covari- 
ant gauge. The gluon propagator can be written 
as 

^ = ~(^ + (l-A)^ 
so that, for A = 1, we recover the Feynman 

gauge. In the general case the presence of a term k^ky leads to an extra term in Eq. (A.3) 
of the form ft(k ■ p\ — k ■ p 2 )- It can easily be seen that this extra term vanishes if one 
replaces 

^ = qA+qh (A.6) 
and makes use of the Dirac equation for massless quarks, 





(A.5) 



Figure 13: Diagram of squark pair produc- 
tion via a s-channel gluon and their subse- 
quent decay. 



qh u(q 2 ) = 0, v( qi ) qA=0. (A.7) 

Next, we turn to the case of non-covariant gauges. We consider the following structure for 
the gluon propagator: 

iT'/j.ki/ + n u k b "~ ? 



D 



1 



n ■ k 



i'"fi ti k^k v 



(A.t 



(n • k) 2 1 

where n can be a time-like, space-like or light-like vector. 

The third term of Eq. (A. 8) vanishes with the same argument as for covariant gauges, as 
well as the term ~ riyk^. 
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The remaining term can be written as 

fi(k-pi-k- p 2 ) = fi(qi ■ Pi + qi ■ Pi - qi ■ Pi - qi ■ Pi) ■ (A. 9) 

Momentum conservation in the on-shell limit implies 

qi-pi = Q2-P2, q2-Pi = qi-P2, (A. 10) 

and therefore the additional factor in Eq. (A. 9) is zero. 
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